On the renormalization of periodic potentials 

I. Nandori^, J. Polonyi^^'^ and K. Sailer^ 

^ Department for Theoretical Physics, Kossuth Lajos University, Debrecen, Hungary 

^ Institute for Theoretical Physics, Louis Pasteur University, Strasbourg, France 

^ Department of Atomic Physics, Lorand Eotvos University, Budapest, Hungary 

(February 1, 2008) 

The renormalization of the periodic potential is investigated in the framework of the Euclidean 
one-component scalar field theory by means of the differential RG approach. Some known results 
about the sine-Gordon model are recovered in an extremely simple manner. There are two phases, an 
ordered one with asymptotical freedom and a disordered one where the model is non-renormalizable 
and trivial. The order parameter of the periodicity, the winding number, indicates spontaneous 
symmetry breaking in the ordered phase where the fundamental group symmetry is broken and the 
solitons acquire dynamical stability. It is argued that the periodicity and the convexity are so strong 
constraints on the effective potential that it always becomes flat. This flattening is reproduced by 
?C ' integrating out the RG equation. 
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^/l • I. INTRODUCTION 

T-H ■ 

^'f-^ I The dynamics generated by a periodic potential represents a challenge in quantum field theory where the usual 

^ i strategies to obtain the solution are based on the Taylor expansion which violate the essential symmetry of the 

t~^ ' problem, the periodicity. One expects two kinds of problems. A local one, the perturbation expansion for the small 

^«0 fluctuations around a minima of the potential should deal with infinitely many vertices in order to preserve the 

' ' periodicity. A global problem appears in the construction of the effective potential as an instability, a conflict between 

^"^ ' the periodicity and the convexity which can be resolved in a trivial manner only, being the constant the only function 

^^ ■ which is periodic and convex in the same time. Another global problem, made rather puzzling by the flatness of the 

^v [ effective potential, is to find an order parameter distinguishing the phase with periodicity from the one where the 

"^ i periodicity is broken spontaneously. 

'"pi ' The goal of this paper is to give a brief presentation of these issues in the case of a two dimensional scalar model. 

JL The pertinent features of the sine-Gordon model, the simplest realization of a periodicity, and the more detailed 

ij ' goals of this paper are presented in section 2. Section 3 contains a very brief introduction into the Wegner-Houghton 

(— I I equation, applied for the periodic potential in section 4. The linearized renormalization group flow is given around 

1" the two-dimensional Gaussian fixed point and the periodic operators of the potential are classified in section 5. The 

. !^ , nonlinear flow is discussed in section 6 by means of numerical integration of the renormalization group equation. 

Si^ ■ Section 7 contains a few remarks about the signatures of the spontaneous breakdown of the fundamental group 

%-{ ' symmetry. Finally, section 8 is for the summary of our flndings. 

a ' 

_ _ I 

II. THE SINE-GORDON, X-Y AND FERMIONIC MODELS 

The simplest example for periodic potential is the two-dimensional sine-Gordon model which is described by the 
Lagrangian 

L = ^id,ci^r+vi^), (1) 

with 

V{(t))=ucosP4>. (2) 

The dynamics and the renormalization have been discussed by means of the straight perturbation expansion and the 
mappings between the sine-Gordon, the Thirring and the X-Y planar models |]l|- [^. The duality transformation 
connecting the sine-Gordon and the X-Y models provides the renormalization group flow in the complete coupling 
constant space [p|jlCt|. 



The mappings between the different models are not exact. In fact, the bosonization of the Thirring model ||l[ and 
the Coulomb-gas representation of a periodic potential |Q are obtained perturbatively. The duality is established 
up to irrelevant terms in the action iQ. But an exact equivalence exists between the X-Y and the compactified 
sine-Gordon model which is obtained by expressing (|l|) in terms of the compact variable |11| 

z{x) = e^^("). (3) 

Such a parameterization makes the kinetic energy periodic, 



l^d,z*d,z + ^{ 



L^^d,z*d,z + -{z + z*) (4) 



and introduces vortices in the dynamics. The models (n|) and (0) are equivalent in any order of the perturbation expan- 
sion in continuous space-time. The modification of the regulator amounts to the introduction of non-renormalizable 
terms in the action which are irrelevant in the UV scaling regime. Thus the equivalence is reached asymptotically 
only when the cutoff is removed. In lattice regularization the theory (|j) is in the same universality class as the X-Y 
model which is described by the action 

s = ^ Y. cos{e^~e.,,) + Y,hcos{e^), (5) 

with T = 0^ . The renormalization of the X-Y model induces a non-trivial value for the vortex fugacity y which 
appears as an additional evolving coupling constant in the compactified sine-Gordon model ||ll|. Therefore the 
complete coupling constant space contains three coupling constants, the external field /i, the temperature T, and the 
vortex fugacity y. 

As it is well known, there are two phases in the y — T plane connected by the Kosterlitz-Thouless transition [^. 
In the low temperature, molecular phase T < Tkt the vortices and anti-vortices form closely bound pairs while 
above the transition temperature T > Tkti in the ionised phase they dissociate into a plasma. Due to the duality 
transformation the corresponding two phases appear in the h — T plane as well as those of the dual electric Coulomb 
gas (DCG). In the h — T plane the ionised phase of DCG is realised at low temperature, i.e. at weak coupling 0^ < 0^ 
whereas the molecular phase of DCG is positioned at high temperature, i.e. at strong coupling 0^ > 01- Instead of 
this duality transformation connecting the X-Y and the sine-Gordon models we shall rely in this paper on the more 
direct equivalence of ^ and (||) as noted above [|ll|, which identifies the weakly coupled (small j3) and the strongly 
coupled (large (3) phase of the sine-Gordon model with the molecular and the ionized phase of the X-Y model. The 
vortices of (0) and (0) correspond to each other in this scheme. 

The transition between the phases can be characterized either perturbatively or by means of the vortex dynamics. 
As far as the perturbation expansion is concerned, the normal ordering is sufficient to remove the UV divergences 
for weak couplings. Though the divergence structure of the individual graphs is the same in either phase the partial 
resummation of the perturbation expansion in u produces new UV divergences for 0^ > Sir ||l|,yj^. A further UV 
divergence was found at /3^ = 47r |5|. The double expansion in u and f3'^ /Sn — 1 0] indicates no special singularities 
at /3^ — An and shows that the adjustment of u and the introduction of a wave function renormalization constant for 
the field (f) in the strong coupling phase is sufficient to remove the UV divergences. 

The inspection of the vortex dynamics allows us to follow the transition line for larger values of u. One way the 
vortices arise is that the vertices appearing in the perturbation expansion in (0) form a gas of vortices described by the 
X-Y model 0. Another way to identify the vortices is to use the equivalence of (^ and ^ in lattice regularization 
when the continuum limit is approached. Both ways indicate that the weak and the strong coupling phases are the 
molecular or the ionised phases, respectively from the point of view of the vortex gas and that the phases are separated 
by the Kosterlitz-Thouless transition. It is worthwhile noting that this transition is rather peculiar because (i) its 
driving force, the vortex dynamics is generated by the UV modes rather than the IR ones as in the case of spontaneous 
symmetry breaking, and (ii) it is a higher than second order since the correlation length is infinite in the molecular 
phase and diverges faster than any power of the reduced temperature (exponentially) , therefore one can not introduce 
the critical exponent v in the usual manner. 

The sine-Gordon model possesses a topological current, 

jf^{x) = —ef,^d,,(t){x) (6) 

which is conserved in the semi-classical expansion, when the path integral is saturated by field configurations with 
analytic space-time dependence. The flux defined by j^ is the vortex number, called the vorticity and the soliton 



number in the same time. In this manner the world hnes of the sine-Gordon soHtons end at the X-Y model vortices, 
making the soliton unstable and the topological current anomalous in the ionised phase where field configurations 
with singular space-time dependence survive the removal of the cutoff [|ll|. This destroys or at least modifies the 
bosonization transformation in the high temperature, ionised phase of the vortex gas. In fact, the non-conservation of 
the topological current requires fcrmion number non-conserving terms in the fermionic representation, a fundamental 
violation of the rules inferred from the weak coupling expansion. 

The relation between the sine-Gordon and the X-Y model is summarized in table ||. 

III. DIFFERENTIAL RG APPROACH IN MOMENTUM SPACE 

The challenge in developing an RG method for the sine-Gordon model is that it should follow the mixing of all the 
operators which become relevant in either phase. Since there are infinitely many relevant o per ators in two dimensions 
one needs the functional form of the evolution equations for the blocked Wilsonian action ||lj]- Jl^. We shall use the 
leading order gradient expansion in the Wegner-Houghton equation [ [L2[ to study the renormalization group flow of 
a generalized model with an arbitrary periodic potential. Such a drastic truncation leaves some doubts due to the 
supposed role of the wavefunction renormalization constant in the ionised phase. But the higher order contribution of 
the gradient expansion can only be treated consistently by the use of the effective action instead of the bare one. We 
shall find that the instability arising from the periodicity of the effective potential makes the Legendre transformation 
highly nontrivial and prevents us to use the effective action in the infrared regime of the molecular phase. 

The differential RG transformations are realised by integrating out the high-frequency Fourier components of the 
field variable, in infinitesimal steps in momentum space successively from the UV cut-off k to k — Sk, 



'Sk-Sk[' 



Dl 



^Sfc[0+0'] 



(7) 



where Sk[4>\ stands for the blocked action with cutoff k and the field variables 4> a-nd </)' contain Fourier components 
with momenta p < k — 5k, and k — 6k < p < k, respectively. In every infinitesimal step, the path integration in 
Eq.(|7|) is evaluated by the help of the saddle point approximation. If the saddle point is at cj)' = 0, one can find an 
integro-differential equation for the blocked action, called the Wegner-Houghton equation jlj] . In the local potential 
approximation one uses the leading order expression for the action in the gradient expansion, 



Sk= d^' 



\{d,<l>{x)f 



VkWx)) 



(8) 



and the Wegner-Houghton equation reduces to a differential equation for the scale dependent potential Vk{4>) [p6| 

kdkVk{^) = -fc'^aln(fc2 +a^l4(0)) , (9) 



with a — ir2d(27r)~ , and the solid angle Vtd in dimension d. 

Notice that the argument of the logarithm in ffl) must be non-negative for the expansion made around a stable 
saddle point. If the argument becomes negative at a critical value k^r > 0, given by /c^,. = —d'i,Vk^^{<j>) then the 
Wegner-Houghton equation looses its validity for k < kcr and the saddle point becomes non-zero and the tree-level 
blocking relation must be used. We simplify the saddle point structure of the blocking (m) by retaining the plane 
waves only and we find the evolution equation 



Vk-Sk{4>) = niin 



*V + i 



where p is the amplitude of the plane wave [pi| . 



duVk{4> + 2pcos(7rw)) 



(10) 



IV. PERIODICITY OF THE POTENTIAL 



The symmetry of the action under the transformation 

(j){x) -> 0(x) + A 



(11) 



is to be preserved by the blocking and the potential 14(<I>) must be periodic with period length A. It is actually 
obvious that the blocking, the transformation 



kVk-skiq^) = kVki^) + [k^aln [k" + 52^^(0))] 5k (12) 

preserves that periodicity of the potential. The impact of the periodicity of the potential on the dynamics can be 
understood by recalling that T4($) tends to the effective potential, Ve//(<&) in the IR limit, fc ^ 0. This is because 
both potential give the action density within the functional space $ — {(j){x)). The Legendre transform imposes 
the condition of the convexity on the effective potential [|l9[ and thereby on Vk=o{^)- Since the only periodic and 
convex function is the constant there is no way to have nontrivial effective potential when the transformation (O) is 
a formal symmetry of the action. Notice that this statement holds for any dimensions. The detailed RG study of the 
sine-Gordon model is pursued here to check such a general conjecture in a simpler case. 

One could object by pointing out that it is enough to implement the horizontal Maxwell cut of the effective potential 
between the inflexion points V^'fA^) = 0, where the function is concave. But the problem left by this construction 
is just at the inflection point where the potential obtained in this manner has singular higher derivatives. This is 
unacceptable, as one can see by placing the system in a large but finite box. On one hand, the effective potential 
should be close enough to the one obtained in the thermodynamical limit, i.e. the Maxwell cut should be present up to 
finite size corrections. On the other hand, the effective potential should be regular in the absence of IR divergences, ft 
is easy to see that the singularity at the inflexion point gives rise to smooth but concave region for finite systems after 
being rounded off and the Maxwell cut must be extended further. In the case of the 0^ model where the symmetry 
which is broken spontaneously isj0 -^ —(j) such an extension of the Maxwell cut leads to the degeneracy of the effective 
potential between the minima |2^. When the symmetry (|l^) is broken then infinitely many minima are expected 
in the effective potential. The only Maxwell cut which is stable in the thermodynamical limit is thus the one which 
renders Veff{^) constant. 

In order to allow a convenient truncation of the potential for the numerical solution of the evolution equation and 
to preserve the periodicity we write the Vk{(j)) as a Fourier series. 



oo 



Ffc(0) = ^u„(fc)cos(n/3</)). (13) 

n=0 

For the sake of simplicity we consider only potentials with Z(2) symmetry, T4(0) — Vfc(— 0). The whole scale 
dependence occurs in the Fourier amplitudes Un{k), the 'coupling constants' of the scale dependent potential. In the 
case of a non-trivial saddle point Eq. (|lO| ) can be rewritten as 



Vfc-5fc(0) = min 



p 



oc 

k^P^ + X! ^"(^) cos{nf3(j))jQ{2nf3p) 



(14) 



where Jq stands for the Bessel function and (3 = 27r/A. Let Pkifp) denote the position of the minimum of the bracket 
on the r.h.s. of Eq.(|l4|). Then pk{(t>)js periodic pk{4' + A) = pk{(t>)i since one has to minimise the same expression of 
p for (j) and </> -I- A. Note that Eq. (O) preserves the periodicity. 

The retaining of the higher order contributions in the gradient expansion, among them the leading order being the 
inclusion of a wavefunction renormalization constant Zk{(l>{x)) into the kinetic energy of the action (q) changes the 
situation. The action keeps its period length in the bare field under any circumstance. But the period length in terms 

of the renormalized field (j)R.k{x) = Z ^ {if)Q)4>{x) , where 4>o minimizes the potential Vfc(0), is AZ^((/)o). Such sub- 
leading contributions of the gradient expansion are neglected in the present work. We believe that their contribution 
will not change our results qualitatively. 

It is easier to use the derivative of Eq. (g) with respect to 4> rather than the original equation itself. For dimension 
d the general form of the evolution equation reads as follows, 

1 ^ 

a/32fc'^-2^V(fc) - (d + fcafc)«„(fc)- -/3A:'*-2'^A„p(fc)(d + A:9fe)vp(A:), (15) 

where N is the truncation in the Fourier series, Vn{k) — n(3un{k) and 

Anp{k) = (n - p)vn^pQ{n > p) (16) 

+ (p - n)wp_„6(p > n) 

- {n+ p)wp+„9(iV >n+p), 

with 9(n > n') = {1 if n > n', if n < n'}. 

Reaching the critical value kcr one has to change automatically from the system of equations dl3) to Eq. (n4) . In 
every step Sk in the momentum space, the potential at the scale k — 6k is then found by minimising the expression on 
the r.h.s. of Eq.(nj). After the minimisation, the potential Vk-Sk{4') is expanded in Fourier-series to define the new 
Fourier amplitudes at the scale k — Sk. One repeats this algorithm step by step till fc = 0. 



V. LINEARIZED SOLUTION 

According to the power counting, theories with polynomial interactions are super-renormalizable in dimension 
d = 2. Furthermore, the super-renormalizable interactions correspond to relevant operators in the UV scaling regime. 
How can we have new UV divergences in the ionised phase when the set of the renormalizable operators is fixed? 
The source of the complication is that in the usual perturbative proof of the renormalizability each monomial vertex 
is treated independently. This strategy is sufficient for polynomial interactions but is not necessarily applicable for 
periodic potentials where the symmetry is destroyed by any truncation of the Taylor expansion. 

The treatment of an infinite series of operators instead of a single monomial may cause complications, an impression 
of having renormalizcd a manifestly non-renormalizable model. In fact, we may find an infinite series of irrelevant 
operators in a renormalizable model, showing the possibility of the removal of the cutoff in the presence of non- 
renormalizable operators. An obvious example is when a regulator, represented as an interaction vertex, yields 
irrelevant operators. We find this situation by introducing the finite difference operator appearing in the lattice 
regularized theories in the continuum. The finite difference operator generates an infinite power series of the gradient 
whose monomials are non-renormalizable. The basic question in the removal of the cutoff is whether the series of 
the irrelevant operators is chosen in such a manner that the divergences can be removed by the fine tuning of a 
finite number of parameters in the action. The infinite series of irrelevant operators is usually required by some 
symmetry of the theory, such as the periodicity in momentum space on the lattice |24|, the global 0{N) symmetry of 
the nonlinear sigma model |2^ or gauge theories on the lattice |2^ . The symmetry imposes such constraints on the 
radiative corrections that the divergences can in fact be removed within the given functional family of the action and 
the apparently non-renormalizable model becomes renormalizable. 

One touches upon here a fundamental difference between the ways renormalization group is used in statistical 
mechanics and particle physics. In statistical mechanics the UV cutoff is physical and we may not ignore the effects 
taking place at that scale. In particle physics we insist that the cutoff is sufficiently far from the scale of the 
phenomenon we are interested even in the effective theories. The gain coming from this constraining of the set of 
observables is that universality arguments apply and it is enough to consider renormalizable modelaj. But we may 
need irrelevant operators in statistical mechanics which lead to complicated non-renormalizable models. The error 
caused by the omission of the irrelevant pieces in the particle physical applications is negligible. We can turn this 
insensitivity into a freedom: The suppression mechanism responsible of this simplification gives us the possibility to 
include irrelevant pieces in the theories without specific fine tuning of their coupling constants in particle physics so 
long as the UV divergences can be removed. 

In the case of the periodic potential we have the opposite effect, a restriction on the renormalizability due to the 
presence of infinitely many vertices in the model. This is because the new UV divergences of the ionised phase arise 
from the summation of infinitely many, individually finite graphs B, B , ||]. We show that this highly nontrivial 
effect can be reproduced in a very simple manner in the framework of the functional form of the renormalization group 
method. We shall use the Fourier amplitudes u„(fc) as coupling constants and we linearize the d = 2 dimensional 



renormalization group flow ( |15[ ) around the fixed point u„ = by assuming \d^Vk{(j))\ ^ k^- Notice that the UV 
Gaussian fixed point is well-defined in any dimensions because modes with non- vanishing momentum are considered in 
the blocking. The singularity building up at or below the lower critical dimension can be found in the flow generated by 
such a blocking procedure as an instability and as an inconsistency of the massless IR fixed point without influencing 
the UV scaling laws. 

It is more natural to express the flow in terms of the dimensionless coupling constants, u{k) — k'~'^u{k). The 
solution of the linearized renormalization group equation satisfying the initial conditions M„(fc = A) = u„a is 



Un{k) = UnA i-TJ ■ (17) 

Thus, the coupling constants Un{k) are relevant, marginal, or irrelevant for /J^ < Sn/n'^, 0^ — Stt/ti^, or 0^ > Stt/ti^, 
respectively. 

For /J^ > Stt all dimensionless coupling constants are irrelevant and ( p7| ) is consistent and keeps the trivial saddle 
point of the blocking stable. This result indicates the inaccessibility of the Gaussian fixed point in the UV limit, the 
non-renormalizability of the model. The infrared fixed point is a trivial, non-interacting massless theory. 



^The universality might need certain generalization in case of several scaling regimes or instabilities 



It is instructive to compare this result with the 0{u^) prediction of the double expansion in u and /3^/87r — 1 
of the sine-Gordon model which shows the possibility of absorbing the UV divergences into u and a wavefuntion 
renormalization constant for the field ^ . The doubt that the higher order contributions of the perturbation expansion 
may render the theory non-renormalizable is resolved in ref. g by relying on the bosonization method. Unfortunately 
there are problems with the bosonization in the ionised phase which is not surprising since this mapping is based 
either on the saddle point or the weak coupling expansion. The problem is the anomaly of the topological current 
which requires the presence of fcrmionic operators with odd power in the fermion field pT| . We believe that the theory 
is non-renormalizable in this phase but the real proof requires further efforts. 

The potential becomes flat under RG transformation in the limit A: — > and the flat effective potential suggests the 
presence of a massless particle which is in contradiction with the dimensionality of the system ||23| . The resolution of 
this apparent contradiction is that according to (O) the non-Gaussian vertices of the model tend to zero sufhcicntly 
fast in the infrared limit to suppress the IR divergences of the perturbation expansion when each mode in the loop 
integral is coupled to the effective coupling constants at the appropriate scale. 

The critical point at 0^ = Stt is a well-known result for the sine-Gordon model. Furthermore it was known that the 
higher harmonics, corresponding to the vortices with higher vorticity are irrelevant around the critical point [^. What 
is interesting in our solution is that one sees a change in the scaling laws for the higher harmonics at a finite distance 
away from this critical point. In fact, the n-th harmonic is found to be relevant for 0^ < Sir/n^. The apparency 
of a new relevant operator implies not only a new renormalization condition but the possibility of new ultraviolet 
divergences. Though we find no evidence for the singularity at /3^ = 47r [g, a series of new singularities is expected 
at /3^ = Sir/n^, due to the vortices with higher vorticity. 

VI. NUMERICAL SOLUTION 

For /3^ < Stt the first few Fourier-amplitudes are relevant, that is they increase for decreasing value of k and 
consequently the linearisation ceases to be reliable. The solution can only be found numerically in this case. Monitoring 
the coupling constants numerically we compare the solution of Eq. ( |l5| ) with the result obtained from the analogous 
relation for the polynomial potential |2^ , 

Vk{q^)^y^^gnr- (18) 

The initial conditions for the polynomial potential were chosen g2 — —0.001, 54 — 0.01 and .g„ = 0, if n > 4. 
Then the saddle point remains trivial for any k and the singularity of Eq.(||) at fc^^(0) — —d'^Vk{(j)) is avoided. In 
order to compare the periodic and the polynomial cases we choose the initial conditions for the Fourier-amplitudes 
Unik = A) for /3^ = 0.1/3^ such, that after Taylor-expansion the initial conditions for the polynomial case are recovered. 
Therefore, the initial conditions for various truncations N of the Fourier-series are different. For the increasing values 
of N the coupling constants gn{k) determined from the periodic potential by the Taylor-expansion approach the 
running coupling constants of the polynomial potential in the UV regime. This is understandable since the quantum 
field (j) does not feel the global properties of the potential, i.e. the periodicity due to its small fluctuations in the UV 
regime. Just the opposite holds in the IR regime, where the field fluctuations become larger and they make the global 
features of the potential manifest. Therefore, it is expected that the solutions for the periodic and the polynomial 
potentials become different in the IR regime. 

We integrated numerically Eq.(p^) starting from the UV cutoff A = 1 down to the critical value k1^{(j}) = —dV/k^^ {4>) 
by using the fourth-order Runge-Kutta method and 5k — lO^^fc with p = i oy p — A. There were no changes in 
the numerical results by increasing p further. In FigJ^ and Fig.0 we show the scaling of the dimensionful coupling 
constants (72 and (74 for different truncations N . Increasing the value of N the differences in the results obtained for 
the periodic potential decrease. 

There are relevant coupling constants for 0^ < Stt which become large enough to destabilize the trivial saddle point 
of the blocking and we reach a non- vanishing critical value k^r where the saddle point becomes non-trivial and the 
tree-level blocking Eq.(|lO|) must be used. By following the solution of this equation all dimensionful Fourier-amplitudes 
are found to approach zero as fc ^ 0. The typical behavior is depicted in Fig.|] and Fig.g. 

There is a remarkable difference in the behavior of the theory with a periodic potential and that with the cor- 
responding polynomial potential p9]. Namely, that all the dimensionful coupling constants g„(fc) obtained for the 
periodic potential tend to zero in contrary to those of the polynomial potential which remain finite as fc ^ 0. In Fig.|3| 
we show this flattening startin g fr om the value kcr- 

Integrating numerically Eq.dlj) we have shown that under the successive infinitesimal RG transformations the 
periodic potential becomes a constant potential up to the accuracy 10~^. We compared our results with the analytic 
formula for the saddle point amplitude Pk{4>) obtained for the polynomial potential El: 



Pfc('/') = ^(0™c(fc)-|0|). (19) 

For periodic potentials V{^ + A) = V7$), the amplitude pk{<P) should be periodic in the field variable with the length 



of period A for any scale k (see Eq.(14)). Thus, for periodic potentials, an expression similar to Eq.(19) is valid in 
the period G [—A/2, A/2] and then the same pattern of the function pk{4') is repeated in all other periods. In the 
particular case investigated by us the minus sign in expression (n9[) should be changed to a plus sign and the field 
independent term is 4>?,Zci^c7-) = —A/2 and (j^^ZA^) — 0. Therefore we can compare the result p^ZoW — ^1^1 ^'^^^ 
that obtained by numerical integration of Eq. (1#, see Fig.0. We have established that with the increasing number 
of Fourier-modes taken into account the computed curves get closer to the dashed line defined by p'^'ZoW — il^l- 

VII. SPONTANEOUS SYMMETRY BREAKING 

We found two different phases of the model, /3^ < Sir and /3^ > Stt, for small u with different scaling laws. The 
effective potential flattens out in either case. How can we reconcile this result with the usual perturbation expansion 
of the sine-Gordon model ||l| where one expects small fluctuations around one of the minima of the periodic potential 
for /3^ < Stt? What does happen with the symmetry ^ (/) -I- 27r? The flattening of the effective potential makes the 
naive order parameter to detect the realization of the periodicity, / dx(j){x), useless. What is the corresponding order 
parameter? 

The common difficulty behind these questions is that the physical configuration space is actually multiply connected. 
In fact, the sine-Gordon Hamiltonian describes a family of coupled pendulums and the configurations (j){x) and (/)(x)-|-A 
are physically indistinguishable. We may use either the multiply connected description in terms of the complex variable 
g27ri0(a;)/A ^ [/(I) Or the Unconstrained (j){x) g TZ variable on the covering space for the description of the system, and 
the transformation ( pT| ) belongs to the fundamental group tti{U{1)) = Z. 

The special difficulty in detecting the spontaneous breakdown of the fundamental group symmetry is the degeneracy 
of the vacuum energy density, the only resolution of the conflict between the periodicity and the convexity of the 
effective potential in either phase. The flattening of the effective action in the region of the Maxwell-cut indicates the 
weakness of the restoring force acting on the fluctuations around the equilibrium position and raises the possibility 
of the collapse of the topological stability of the solitons. This, in turn, suggests the winding number as an order 
parameter to distinguish between the phases with the usual and the unusual realization of the fundamental group 
symmetry. This circumstance reflects a further similarity with gauge models at finite temperature ||22| . 

To understand these questions better we start by distinguishing the local potential of the effective theory Vk{(j)) 
with a low cutoff k from the effective potential, Veff{4>) = Vfc=o ('/')• In theories with infrared stable dynamics the limit 
k —> is safe and this difference is negligible. But the infrared instability inducing the mixed phase and the Maxwell 
cut in the effective potential makes the limit A: — > more involved even for theories with massive particles only pi[ . 
By turning this complication into an advantage, we suggest to consider the infrared instabilities as the signature of 
the spontaneous symmetry breaking. 

Let us consider the dynamics of the modes with momentum p < k ^ 0. The local potential Vk{4>) of the corre- 
sponding effective theory is non-trivial and periodic and the perturbation expansion might be justified around one of 
the minima for k > 0. The real question is the relative speed the different coupling constants approach as /c ^ 0, 
and whether the saddle point and the perturbation expansion remain consistent in this limit. As mentioned above, 
we adopt the infrared instability, i.e. the stability of the trivial saddle point of the blocking as a signature of the 
absence of the change of the symmetry pattern of the model. 

The saddle point of the blocking remains trivial in the ionic phase and we expect no spontaneous symmetry 
breaking there. The coupling constants decrease as fc — > 0, the theory becomes trivial as it happens with other 
non-renormalizable models. The infinitesimal fluctuations gradually flU up the valleys of the local potential and the 
naive order parameter for the periodicity, J dx4>{x), decouples. 

The large amplitude, tree level fluctuations which lie beyond the realm of the perturbation expansion but are picked 
up by the tree-level renormalization start to fill up the valleys of the potential in the molecular phase already at a 
finite scale, k < kcr- We take this instability as the indication of spontaneous symmetry breaking. This seems to be 
in agreement with the perturbative approach which is meaningful only when the potential has a non-trivial structure. 

In order to identify the order parameter for the periodicity we elucidate the topological differences between the 
two phases. Let us introduce periodic boundary conditions in the Euclidean space-time and consider the following 
quantities. One of them is the winding number, 

Q(xO)= [dx'jo{x°,x^), (20) 



the space integral of the topological current density at a given time, 1°. We shall argue that this is the order parameter 
for the fundamental group symmetry. One can construct another, independent topological invariant by exchanging 
the time and the space axes. But this does not modify the discussion what follows. The other quantities are the 
vorticities of the space-time regions before and after the time x^ , 

V±= f d?vd^j^{y). (21) 

J±(x°~y°)>0 

The flux of the topological current agrees with the vorticity of the enclosed vortices Q = V- = —V+ |ll| ]. 

In the molecular phase the distance between the vortex-anti vortex pairs is shrinking with the increasing UV 
cutoff and the quantum fluctuations cannot change the value of Q. This makes the path integration consistent 
when constrained within a homotopy class, characterized by a fixed value of Q. In fact, the consistency of the path 
integration is the requirement that the path integral as the function of the end point of the trajectories should satisfy 
the (functional) Schrodinger equation. The Schrodinger equation can be derived for the path integral by performing 
infinitesimal variations on the trajectories at the final time. Thus the path integral is consistent if the functional space 
over which it is evaluated is closed with respect to infinitesimal deformation of the trajectories. Since the discontinuous 
field configurations are suppressed in the path integral the quantization process is well defined in a given homotopy 
class [2^]. Note that the consistent constraining the functional integration into a given homotopy class removes the 
fundamental group symmetry. This is obvious in the semi-classical quantization of a soliton where the dynamical 
stability of the whole construction comes from the spontaneous breakdown, the suppression of the winding number 
changing processes. Thus the condition for the stability of the solitons, the sufficient smoothness of the configurations 
dominating the path integral, is in the same time the signature of the breakdown of the fundamental group symmetry. 

There are vortex-anti vortex pairs with cutoff independent separation in the ionised phase and Q fluctuates in an 
uncontrollable manner. The path integration cannot be constrained into a given homotopy sector, the fundamental 
group symmetry is realized but the semi-classical structure based on smoothness is destroyed, the solitons become 
unstable. In the same time the bosonization relations are either lost or fundamentally modified because the soliton 
(fermion number) non-conserving processes require the introduction of operators with odd powers of the fermion field 
in the action. 

The susceptibility of the topological charge, 

X = (0') - {Q? (22) 

may serve as a disorder parameter to distinguish the different realizations of the periodicity. In fact, the vortex 
fugacity tends to zero or stays finite in the molecular phase or the ionized phase when the cutoff is removed according 
to the Kosterlitz-Thouless RG flow. Since the topological susceptibility is vanishing whenever the fugacity is zero, 
the former may serve as an order parameter. 

VIII. SUMMARY 

A simple two dimensional scalar model with periodic potential was investigated in this paper. The first question 
considered was the renormalization of the potential. The study of this problem requires the handling of infinitely 
many operators which was achieved by the Wegner-Houghton equation. We found a disordered phase where the model 
is non-renormalizable and trivial. In another, ordered phase the relevant operators compatible with the symmetry 
were identified. 

The periodicity and the convexity impose triviality on the effective potential, a phenomenon verified in detail in 
both phases. The coupling constants approach zero regularly in the disordered phase as the theory becomes trivial, 
leading to the flattening of the effective potential. Instabilities were found in the ordered phase where the effective 
potential is flattened out by the Maxwell construction only, indicating a nontrivial dynamics behind the trivial final 
result. 

It was pointed out that the configuration space of the model is multiply connected and that the winding number 
is suggested as a non-local order parameter distinguishing the explicit realization and the spontaneous breakdown of 
the fundamental group, the shift of the field variable by the period length of the potential. The relation between the 
applicability of the semi-classical arguments and the spontaneous breakdown of the fundamental group symmetry is 
shown in the context of our model. 

We have started to extend this work to include the wavefunction renormalization constant for the field, and to 
investigate the phase structure in higher dimensions. The scalar model with periodic potential can be considered 
as a non-linear 0(2) model with a symmetry breaking term, and one expects similarities with non-Abelian gauge 



theories due to the compact nature of the dynamical variable. In fact, the molecular phase is asymptotically free and 
the vacuum is filled up with vortices in its large distance structure. Furthermore, the spontaneous breakdown of the 
fundamental group symmetry is reminiscent of the breakdown of the center symmetry in finite temperature gauge 
theories. Thus, our results offer interesting lessons to be learned in constructing the quark confinement mechanism. 

We ignored the wavefunction rcnormalization constant in this work. In order to go beyond this approximation 
and to take into account some higher order contributions in the gradient expansion one has to turn to the evolution 
equations for the effective action, instead of the bare action ||l5[, ||l^, ||l^. But the problem is just that the effective 
action, appearing in the rcnormalization group treatment, hides a large part of the dynamics due to the Maxwell cut. 
This is because the effective action, obtained by a smooth cutoff is after all the Legendre transform of the logarithm of 
a path integral which does contain the IR modes. These modes, though they appear with small amplitude, generate 
the Maxwell cut. It is not clear to us how to improve the renormalization group method to make it applicable for 
models with spontaneous symmetry breaking or with compact variables (e.g. gauge models with gauge fixing which 
are based on compact gauge group) beyond the local potential approximation. 
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FIG. 1. Comparison of the scaling of the dimensionful coupling constant g2 obtained for the polynomial (dashed-line with 
stars) and for the periodic potential for various values of the truncation A^ = 2, 3, 10, 20, 30 in the Fourier series. The figures on 
the left (right) show the scaling of the coupling above (below) kcr- Below k^r we show the scaling oi g2 for the cases A'' = 2, 3, 10. 
The increasing order of the full-line curves corresponds to increasing A''. 
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FIG. 2. The scaling of the dimensionful coupling constant 34 in the same cases as for 32 in FigM The figures on the left 
(right) show the scaling of the coupling above (below) kcr- 
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FIG. 3. Flattening of the periodic potential below kcr- The decreasing order of the curves corresponds to decreasing values 
of the scale k with the step Sk = kcr/W- 
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FIG. 4. Comparison of the function pk=o{(i>) obtained by numerical integration of Eq.(|l4|) with the analytic expression 



p'l^Zoi't') — k\4'\ (dashed-line) . 



XY-model with external field 


compactified sine-Gordon model 


external field h 


Fourier amplitude u — h/(3'' 


temperature T 


coupling constant /3 


molecular phase of the vortex gas T <Tc 


weak coupling phase /3"^ < P^ 


ionised phase of the vortex gas T > Tc 


strong coupling phase /3^ > f3^ 


vortex, anti-vortex 


soliton creation and annihilation 



TABLE I. Comparison of the XY and the sine-Gordon model 
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